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1 Introduction 



Circle patterns as discrete analogs of conformal mappings is a fast developing field of research on the 
border of analysis and geometry. Recent progress in their investigation was initiated by Thurston's 
idea flq] about approximating the Riemann mapping by circle packings. The corresponding con- 



vergence was proven by Rodin and Sullivan [15]. For hexagonal packings, it was established by He 
and Schramm Q that the convergence is C°°. Classical circle packings comprised by disjoint open 
disks were later generalized to circle patterns, where the disks may overlap (see for example M). 
In [p| ; Schramm introduced and investigated circle patterns with the combinatorics of the square 
grid and orthogonal neighboring circles. In particular, a maximum principle for these patterns was 
established which allowed global results to be proven. 

On the other hand, not very much is known about analogs of standard holomorphic functions. 
Doyle constructed a discrete analogue of the exponential map with the hexagonal combinatorics 
||], and the discrete versions of exponential and erf- function, with underlying combinatorics of the 
square grid, were found in |f[q |. The discrete logarithm and z 2 have been conjectured by Schramm 
and Kenyon (see Jl7[). 

In a conformal setting, Schramm's circle patterns are governed by a difference equation which 
turns out to be the stationary Hirota equation (see 16, ||]). This equation is an example of an 



integrable difference equation. It appeared first in a different branch of mathematics - the theory 
of integrable systems (see IS] for a survey). Moreover, it is easy to show that the lattice comprised 



by the centers of the circles of a Schramm's pattern and by their intersection points is a special 
discrete conformal mapping (see Definition p] below). The latter were introduced by || in the 
setting of discrete integrable geometry, originally without any relation to circle patterns. 

The present paper is devoted to the discrete analogue of the function f(z) = z" 1 , suggested first 
in 0. We show that the corresponding Schramm's circle patterns can be naturally described by 
methods developed in the theory of integrable systems. Let us recall the definition of a discrete 
conformal map from ||. 

Definition 1 / : Z 2 — > R 2 = C is called a discrete conformal map if all its elementary quadri- 
laterals are conformal squares, i.e. their cross-ratios are equal to -1: 



{.fn,rmfn+l,mifn+l,m+lifn,m+l) • — 

= -1 (1) 



(fn,m fn+l,m) (/n+l,m+l fn,m+l 



(/n+l,m fn+l ,m+l)(/n ,m+l fn,m) 



This definition is motivated by the following properties: 

1) it is Mobius invariant, and 

2) a smooth map / : D C C — > C is conformal (holomorphic or antiholomorphic) if and only if 
V(x,y)€D 

^Q{f(x,y), f{x + e,y)f(x + e,y + e)f(x,y + e)) = -1. 
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For some examples of discrete conformal maps and for their applications in differential geometry 
of surfaces see || [lC} . 

A naive method to construct a discrete analogue of the function f{z) = z 7 is to start with 
/n,o = nl ■> n > 0, /o jm = (im) 7 , m > 0, and then to compute f n ^ m for any n, m > using equation 
(|l]). But a so determined map has a behavior which is far from that of the usual holomorphic maps. 
Different elementary quadrilaterals overlap (see the left lattice in Fig. [l]). 




Figure 1: Two discrete conformal maps with close initial data n = 0,m = 0. The second lattice 
describes a discrete version of the holomorphic mapping z 2//3 . 



Definition 2 A discrete conformal map f nm i- s called an immersion if the interiors of adjacent 
elementary quadrilaterals (f n ,m, fn+l,m, fn+l ,m-\-l ; fn, 

m+ l) are disjoint. 

To construct an immersed discrete analogue of z 7 , which is the right lattice presented in Fig. |l[ 
a more complicated approach is needed. Equation ([[]) can be supplemented with the following 
nonautonomous constraint: 

/• n (/n+l,m fn,rn)(fn,m fn—l,m) . n (/n,m+l fn,m){fn,m fn,m— l) /n\ 

lfn,m = 2n — — r + 2m — — , (2) 

\Jn+l,m Jn—l,m) \Jn,in+l Jn,m—1) 

which plays a crucial role in this paper. This constraint, as well as its compatibility with (|l]), is 
derived from some monodromy problem (see Section 2). Let us assume < 7 < 2 and denote 
= {(n,m) G Z 2 : n, m > 0}. Motivated by the asymptotics of the constraint (FJ) as?i,m-> 00, 
and by the properties 

z 7 (R+) € R+, z 7 (fR+) e e 7 ^ /2 R+ 
of the holomorphic mapping z 7 , we use the following definition H, ^ of the "discrete" z 7 . 

Definition 3 The discrete conformal map Z 7 : Z 2 — > C, < 7 < 2 is the solution of ^) and 
^) with the initial conditions 

Z 7 (0,0) = 0, Z 7 (1,0) = 1, Z 7 (0,1) = e 7 " /2 . (3) 
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Obviously Z 7 (n,0) £ R + and £7(0, m) G e 77ri / 2 (R+) for any n,m G N. 

Fig. 2 suggests that Z 7 is an immersion. The corresponding theorem is the main result of this 
paper. 

Theorem 1 The discrete map Z 7 for < 7 < 2 is an immersion. 

The proof is based on analysis of geometric and algebraic properties of the corresponding lattices. 




Figure 2: Schramm's circle pattern corresponding to Z 2 / 3 

In Section |H we show that Z 7 corresponds to a circle pattern of Schramm's type. (The circle pattern 
corresponding to Z 2 / 3 is presented in Fig. ||.) Next, analizing the equations for the radii of the 
circles, we show that in order to prove that Z 7 is an immersion it is enough to establish a special 
property of a separatrix solution of the following ordinary difference equation of Painleve type 

/ -.w 2 -1 \ ( x n+l — l%n \ / 2 1 1 \ ( x n~l + \ 

(n + l)(x n - 1) — - n[x n + 1) — = 

\l + X n X n+1 J \l + X n -lX n ) 

Namely, in Section [| it is shown that Z 7 is an immersion if and only if the unitary solution 
x n = e*° n of this equation with xq = e* 77r//4 lies in the sector < a n < ir/2. Similar problems 
have been studied in the setting of the isomonodromic deformation method |11, q. In particular, 
connection formulas were derived. These formulas describe the asymptotics of solutions x n for 
n — ► 00 as a function of xo (see in particular @). These methods seem to be insufficient for our 
purposes since we need to control x n for finite re's as well. The geometric origin of this equation 
permits us to prove the property of the solution x n mentioned above by purely geometric methods. 
Based on results established for Z 7 , we show in Section || how to obtain discrete immersed analogs 
of z 2 and log z as limiting cases of Z 7 with 7^2 and 7 — > 0, respectively. Finally, discrete analogs 
of Z 7 for 7 > 2 are discussed in Section 0. 
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2 Discrete Z 1 via a monodromy problem 

Equation (jlj) is the compatibility condition of the Lax pair 

'I'n+ljm — fra,m^ra,m ^n,m+l — Vn,rn^ n,rn (4) 



found by Nijhoff and Capel [13]: 



U n) m = ( A j Vn,m = ( A ^ j , (5) 

where 

Un,m — fn+l,m fn,mi ^n,m — fn,m+l fn,m- 

Whereas equation (|l]) is invariant with respect to fractional linear transformations f n ^ m — > (pf n ,m + 
q)/( r fn,m + s), the constraint (g) is not. By applying a fractional linear transformation and shifts 
of n and m, (J2J) is generalized to the following form: 

P 2 f , r <-,/ ±\ {fn+l,m fn,m){fn,m fn—l,m) 



0fn,m + lfn,m + 5 = 2(n 



(/n+l,m fn—l,m) 



n r ,\{fn,m+l fn,m){fn,rn fn,m—l) ,„-. 

2{m-ip) 77 —f ^ ' W 

\Jn,m+l Jn,m—l) 

where /3, 7, <5, 0, tp are arbitrary constants. 

Theorem If : Z 2 —* C is a solution to the system (Q, i/ and on/y if i/iere exists a solution 
^n,m to satisfying the following differential equation in X: 

jLxii - A \l> A - Bn ' m I Cn ' m l Dn ' m (7) 

dX 1 + A 1 — A A 

with X— independent matrices B n ,rm Cn,mi D n ^ m . The matrices B n ^ m , C n)m , D n ,m 

in ((J) are of 

the following structure: 

^ / ^n,m %,m"n-l,m \ r 1 j 

u n,m + M n - l,m \ 1 n n—l,m J 2 

_ m-lp ( V n ^ m V„, m W n ,m-l \ _ r 

L'n^ — I -i Jo 

u n,m + ^n,m- 1 V - 1 v n,m—l J £ 

( _2_/3f f2 _7f _i \ 

n — I 4 2-l n ' m 2Jn,m 2-> n > m 2 

n ' m ~ -I 7 , £f 

\ 2 4 ~ r 2-l n > m / 

77ie constraint (j^j is compatible with (j^j. 
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the proof of this theorem is straightforward but involves some computations. It is presented in 

Appendix A. 

Note that the identity 

det tf n , m (A) = (1 + A) n (l - A) m det tf ,o(A) 

for determinants implies 

Tl Tfl 

tr^„, m (A) = — - — + a(A), (8) 

where a(A) is independent of n and m. Thus, up to the term D n ^ m /X, equation (0) is the simplest 
one possible. 

Further, we will deal with the special case in (^) where P = 5 = <j) = ip = 0, leading to the discrete 



Z 7 . The constraint (y) and the corresponding monodromy problem were obtained in [12] for the 
case 7 = 1, and generalized to the case of arbitrary 7 in H. 

3 Circle patterns and Z 1 

In this section we show that Z 1 of Definition is a special case of circle patters with the combina- 



torics of the square grid as defined by Schramm in [16]. 



Lemma 1 Discrete f n , m satisfying and (^) with initial data /o,o = 0, fi$ = 1, /o,i = e * a has 
the equidistant property 

/2n,0 — /2n-l,0 = fin+1,0 _ /2n,0> /o,2m — /o,2m-l = /o,2m+l _ /o,2m 

for any n>l,m>\. 

Proof: For m = or n = the constraint (§) is an ordinary difference equation of the second order. 
The Lemma is proved by induction. 

Given initial /o^, /o,i and fxfi the constraint (||) allows us to compute f n fi and /o, m for all 
n,m> 1. Now using equation (||) one can successively compute fn im for any n,m £ N. Observe 
that if |/n+l,m — fn,m\ = \fn,m+l - fn,m\ then the quadrilateral (f n ,rn,fn+l,m,fn+l 
is of the kite form - it is inscribed in a circle and is symmetric with respect to the diame- 
ter of the circle [f n ,m, fn+i,m+i]- If the angle at the vertex / n>m is tt/2 then the quadrilateral 
(fn,m, fn+i,m, jn+i.m+i, /n,m+i) is of the kite form too. In this case the quadrilateral is symmetric 
with respect to its diagonal [f n ,m+l, fn+l,m]- 

Proposition 1 Let / n . m satisfy and (j^) in with initial data /o,o = 0, /o,i = 1, /o,i = e ia . 
Then all the elementary quadrilaterals (f n ,m, fn+i,m, fn+i,m+i, fn,m+i) are of the kite form. All 
edges at the vertex f n ^ m with n + m = (mod 2) are of the same length 

|/n+l,m fn,m\ — \fn,in+l fn,m\ — \fn— l,m fn,m\ — \ fn,in—l fn,m\- 

All angles between the neighboring edges at the vertex f n)Tn with n + m = 1 (mod 2) are equal to 
tt/2. 
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Proof follows from Lemma |l] and from the above observation by induction. 

Proposition [l] implies that for any n,m such that n + m = (mod 2), the points f n +i,m, fn,m+i, 
fn-i,m, fn,m-i he on a circle with the center / n>TO . 




Figure 3: Discrete conformal maps of Schramm type: sublattices and kite-quadrilaterals, n + m = 
(mod 2) 

Corollary 1 The circumscribed circles of the quadrilaterals (/n-i, mi /n,m-ii /n+i,m, /n,m+i) with 
n + m = (mod 2) form a circle pattern of Schramm type (see fcf(\]), i.e. the circles of neighboring 
quadrilaterals intersect orthogonally and the circles of half-neighboring quadrilaterals with common 
vertex are tangent (see Fig. ^). 

Proof: Consider the sublattice {n, m : n + m = (mod 2)} and denote by V its quadrant 

V = {z = N + iM: N,M £ Z 2 ,M > |JV|}, 

where 

N = (n-m)/2, M=(n + m)/2. 
We will use complex labels z = N + iM for this sublattice. Denote by C(z) the circle of the radius 

— \fn,rn fn+l,m\ — \fn,m fn,m+l\ — \fn,m fn— l,m| — \fn,m fn,m— 1| (9) 

with the center at fN+M,M-N = fn,m- From Proposition [l] it follows that any two circles C(z), 
C(z') with \z — z'\ = 1 intersect orthogonally, and any two circles C(z), C(z') with \z — z'\ = \[2 
are tangent. Thus, the corollary is proved. 

Let {C(z)}, z G V be a circle pattern of Schramm type on the complex plane. Define f n m '■ ~^ C 
as follows: 

a) if n + m = (mod 2) then / n , m is the center of + i^), 
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b) if n + m = 1 (mod 2) then f n>m := + z^±f=±) n C(^f±± + *^±i) = C(^=f±± + 

n C(^=f=^ + i2±f±i). Since all elementary quadrilaterals 
(/n,m,/n+l,m,/n+l,m+l,/n,m+l) are of the kite form, equation @) is satisfied automatically. In 
what follows, the function /n )m , defined as above by a) and b), is called a discrete conformal map 
corresponding to the circle pattern {C(z)} . 

Theorem 3 Let / n . m satisfying |Jj and with initial data /o f o = 0, /o,i = 1, /o,i = e lQ , 6e an 
immersion. Then R(z) defined by f^j satisfies the following equations: 



R{z)R{z + 1)(-2M - 7 ) + + l)R(z + 1 + i)(2(N + 1) - 7 )+ 

+ 1 + + i)(2(M + 1) - 7) + R(z + i)R(z)(-2N - 7) = 0, 

/or z £ Vj := V U {-iV + - 1)|2V € N} and 

(iV + M){R{z) 2 - R(z + l)R(z - i))(R(z + i) + #(2 + 1))+ 

(M - N){R{z) 2 - R{z + i)R(z + l))(R(z + 1) + #0 - i)) = 0, 



(10) 



(11) 



for z G V int := V\{±iV + iJV|/ST G N}. 

Conversely let R{z) : V — > R + satisfy ( fIZ| j /or z 6 V; and ([Tip /or z G Vj n £. T/ien -R(z) defines 
an immersed circle packing with the combinatorics of the square grid. The corresponding discrete 
conformal map f n ^ m is an immersion and satisfies (^). 

Proof: Suppose that the discrete net determined by / n>m is immersed, i.e. the open discs of 
tangent circles do not intersect. Consider n + m = 1 (mod 2) and denote f n +i,m = fn,m + ne 1 ^, 
fn,m+i = fn,m + ^ 2 e i/3 , f n -i, m = fn,m ~ ^2,^ ', /n,m-l = fn,m - ine l/3 , where n > are the radii 
of the corresponding circles. The constraint (|2|) reads as follows 

7/„, m = & (2n^- + 2im^-) . (12) 
V n + r 3 r 2 + r 4 J 

The kite form of elementary quadrilaterals implies 

, _ , _ ip (n - ir 2 f _ ij3 (ri + ir 4 ) 2 

JYi+l/m+1 — Jn+l,m C r\ „ r> , Jrt+1 ,m— 1 — Jn+l.m 6 7*1 9,9 

rf + r 2 rf + r| 



Computing f n +2,m from the constraint (|T^) at the point (n+ 1, m) and inserting it into the identity 
\fn+2,m — fn+i,m\ = n, after some transformations one arrives at 

rir 2 (n + m + 1 — 7) + r 2 r3(— n + m + 1 — 7) + r^r^—n — m + 1 — 7) + r 4 ri(n — m + 1 — 7) =0. (13) 

This equation coincides with (|l0|). 

Now let f n +2,m+l = fn+l,m+l + R\C % ^ , fn+l,m+2 = fn+l,m+l + iR2^ , fn,m+l = fn+l,m+l ~ 

R^ 1 ^ 1 fn+i,m = fn+i,m+i — iRifi ■ Since all elementary quadrilaterals are of the kite form we 
have 



R, = n, R 3 = r 2 , eW = -j e *fa+^ 

rf + r% 



2 
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Substituting these expressions and ( |i~2^ ) into the constraint @ for (n + l,m + 1) and using (fb|), 
we arrive at: 

(n + l)rf (r 2 + r 4 ) + mr 2 (rf - r 2 r 4 ) 



#2 



(n + l)r 2 (r 2 + r 4 ) - m(rf - r 2 r 4 ) 
(m + l)r|(ri + r 3) + nr\{r\ — ^1^*3) 



(m + l)ri(ri + r%) — n{r\ — r\r%) 

These equations together with R4 = n, R3 = r 2 describe the evolution (n, m) — ► (n + 1, m + 1) of 
the crosslike figure formed by / n ,m; /n±l,m> fn,m±i with n + m = 1 (mod 2). The equation for i? 2 
coincides with (0). We have considered internal points 2 G Vj„£, now we consider those that are 
not. Equation @ at z = N + i N and z = -N + i(N - 1), iV G N reads as 

R(±(N + 1) + i(N + 1)) = ^ ^L±T— R(N + iiV). (14) 

The converse claim of the Theorem is based on the following Lemma. 



Lemma 2 Let R(z) : V -» R + safe/y @) /or z£V, and ^ /or z = i¥, M G N. Then R(z) 
satisfies: 

a) equation /or z G V\{iV + iiV|JV G N} ; 
6j equation 



(N + M){R(z) 2 - R(z + i)R(z - l))(R(z - 1) + R(z - »))+ 
(M - N)(R(z) 2 - R(z - l)R(z - i))(R(z + + R(z - 1)) 



(15) 



/or z G V\{-N + iiV|iV G N} 
equation 

2 (i!(iTT) + 7^ + fl^ry + grjgy) g(f + ^ + ~ W ~ g 

^ + 1) + fl(z + i) + - 1) + - ») [ } 

for z G V int . 

Proof of this Lemma is technical and is presented in Appendix B. 

Let R{z) satisfy (|lO],[ll]) then the item c) of Lemma |2| implies that at z G Vj n t equation ([l(]) is 
fulfilled. In JL6| ] it was proven that, given R(z) satisfying flllf), the circle pattern {C(z)} with radii 
of the circles R(z) is immersed. Thus, the discrete conformal map f nim corresponding to {C(z)} is 
an immersion. The item b) of Lemma |2] implies that R(z) satisfies (15) at z = N + iN, JVgN, 
which reads 

R(N - 1 + iN)R(N + i(N + 1)) = R 2 (N + iN). (17) 

This equation implies that the center O of C(N +iN) and two intersection points A, B of C(N+iN) 
with C(N - 1 + iN) and C(N + i(N + 1)) lie on a straight line (see Fig. g). Thus all the points / n>0 
lie on a straight line. Using equation ([H]) at z = N + ^V, one gets by induction that / n>m satisfies 
(|2|) at (n, 0) for any n > 0. Similarly, item a) of Lemma ^, equation (|TT| ) at z = — AT + i/V, A^ G N, 
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Figure 4: Straight line / n ,o 

and equation ( |To| ) at z = —N + i(N — 1), iV G N, imply that f n>m satisfies (H) at (0, to). Now 
Theorem ^ implies that /„ )m satisfies (||) in Z+, and Theorem || is proved. 

Remark. Equation (^) is a discrete analogue of the equation Alog(i?) = in the smooth case. 
Similarly equations and ([H]) can be considered discrete analogs of the equation xR y — yR x = 0, 
and equation (|i~0|) is a discrete analogue of the equation xR x + yR y = (7 — 

From the initial condition (||) we have 

12(0) = 1, i*(i)=tan^p (18) 

Theorem ||| allows us to reformulate the immersion property of the circle lattice completely in terms 
of the system (10, 11). Namely, to prove Theorem |l] one should show that the solution of the system 
(|j~0| , 11) with initial data (18) is positive for all z G V. Equation (|l4|) implies 

iZ( ± JV + ijV) = 7(2 + 7) - (2(JV - 1)+7) (19) 
( + j (2- 7 )(4- 7 )...(2iV- 7 )- (i9j 



Proposition 2 Let i/te solution R(z) of and ^iby in V with initial data 

R(0) = 1, i?00=tan^ 

be positive on the imaginary axis, i.e. R(iM) > for any M S Z+. T/ien -R(z) is positive 
everywhere in V. 

Proof: Since the system of equations for R(z) defined in Theorem || has the symmetry N — > —N, 
it is sufficient to prove the proposition for N > 0. Equation (|l0|) can be rewritten as 



fi(z) J R(z + l)(2M + 7) + i?(z)i?( Z + i)(2iV + 7) 
lZ ' j R(z + l)(2iV + 2 - 7) + R(z + i)(2M + 2 - 7) ' 

For 7 < 2, N > 0, M > 0, and positive R(z), R(z + 1), R(z + i), we get R(z + 1 + i) > 0. Using 
R(N + iN) > for all JVeN, one obtains the conclusion by induction. 
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4 Z 1 and discrete Painleve equation 



Due to Proposition ^ the discrete Z 7 is an immersion if and only if R(iM) > for all MeN, To 
prove the positivity of the radii on the imaginary axis it is more convenient to use equation (|2|) for 
n = m. 



Proposition 3 The map f : —* C satisfying ^) and ^) with initial data /o,o = 0, /o,i 
/o i = e ta is an immersion if and only if the solution x n of the equation 



1. 



(n + 1)04-1) 



nix„ + 1 



(20) 



v l -\- X n X n ^-l J \Z -\- X n —\X n 

with xq = e %a l 2 , is of the form x n = e toln , where a n £ (0,7r/2). 

Proof: Let f n , m be an immersion. Define R n := R(in) > 0, and define a n £ (0,7r/2) through 
fn,n+i ~ fn,n = e 2ja ™ (/ n +i,n - fn,n)- By symmetry, all the points lie on the diagonal arg/„ in = 
a/2. 




Figure 5: Diagonal circles 
Taking into account that all elementary quadrilaterals are of the kite form, one obtains 
/ n+2 , n+1 = e ia / 2 (g n+1 + R n+1 e~ ia ^), f n+ i, n+2 = e ia / 2 (g n+1 + R n+1 e ia ^), 



fn+l,n — e' a ' 2 (ffn+l ~ iR n +l^ 



f n>n+1 = e ia/2 (g n+1 + iRn+ie* 



and 



R n+ i = R n t&na n , (21) 
where g n +i = \ fn+i,n+i\ (see Fig. ||). Now the constraint @ for (n + 1, n + 1) is equivalent to 

g««n _|_ ^gi«„+l \ 



7g n+ i = 2(n + l)R n+1 



_|_ e i(a n +a n+1 ) J 
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Similarly, 

/ gia n _i _|_ igia-n 

jg n = 2nR n 
Putting these expressions into the equality 



_|_ g»(a n -i+a„) I 



I 



9n+l =9n + e ia " (R n + iRn+l) 



and using (21) one obtains (|20| ) with x n = e lOLn . This proves the necessity part. 



Now let us suppose that there is a solution x n = e totn of ( p0| ) with a n S (0, vr/2). This solu- 
tion determines a sequence of orthogonal circles along the diagonal e ia//2 R+, and thus the points 
fn,n, fn±i,n, fn,n±i, for n > 1. Now equation d) determines / n>m in Z^. Since a n £ (0,tt/2), the 
inner parts of the quadrilaterals (/ n ,n> fn+i,n, fn+i,n+i, fn,n+i) on the diagonal, and of the quadri- 
laterals (/ n ,n-i 5 /n+l,n-i) fn+l,n, fn,n) are disjoint. That means that we have positive solution -R(z) 
(rf®,(0) for z = iM, z = 1+iM, N G N. (See the proof of Theorem [|) Given R(iM), equation 
(|10D determines -R(z) for all z S V. Due to Lemma ^, satisfies (^, [TT|) . From Proposition 

Q it follows that R(z) is positive. Theorem |3] implies that the discrete conformal map g n , m corre- 
sponding to the circle pattern {C(z)} determined by R(z) is an immersion and satisfies Since 

9n,n 

= f n . n and = fn±l,n, equation (g) implies f n 

,m — 9n,m- This proves the theorem. 

Remark. Note that although ( p0[ ) is a difference equation of the second order a solution x n of (p0|) 
for n > is determined by its value xq = e ia ^ 2 . From the equation for n = one gets 

x x = X » {XI \\- l \ . (22) 

i(( 7 - 1)^2 + 1) 1 > 



Remark. Equation ( |20| ) is a special case of an equation that has already appeared in the literature, 
although in a completely different context. Namely, it is related to the following discrete Painleve 
equation 

2Cn+l 2£ n 

+ i TT-TT = /i + + Cn+l + Cn+ 



1 — X n+ iX n 1 — X n X n __i 



- ^)(r 2 - 1)X„ + r(l - X2)[I(C„ + Cn+l) + ("l) n (Cn " Cn+1 " 2m)] 



(r + X n )(l +rX f , 



which was considered in [14|, and is called the generalized d-PII equation. The corresponding 

transformation^ is 

U x _ (l + i)(x-i) 

y/2{x + l) 

with Cn = n, r = -\/2, /J, = 0, (( n - Cn+l - 2m) = 0, 7 = (2u - ( n + Cn+l)- 
1 We are thankful to A.Ramani and B.Grammaticos for this identification of the equations. 
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Equation ( |20| ) can be written in the following recurrent form: 

x n +i = f(n,x n -i,x n ) := 

nx~ 2 + i(7 - l ) x n-i x n 1 + (7 - 1) + »(2n + l)^^ + (n + l)z 2 
nx^ - 3(7 - l)x n -ix n + (7 - 1) - i(2n + ^Xn-iXn 1 + (n + l)x« 2 
Obviously, this equation possesses unitary solutions. 

Theorem 4 There exists a unitary solution x n of the equation $2(\ ) with 
x n G Aj\{l,i} E S , Vn > 0, where 

Aj :={jf>\pe[0,n/2]}. 



Proof: Let us study the properties of the function (p(n, x, y) restricted to the torus 
T 2 = S 1 x S 1 = {(x, y):x,yeC, \x\ = \y\ = 1}. 

1. The function (f(n, x, y) is continuous on Aj x Aj Vn > 0. (Continuity on the boundary of Aj x Aj 
is understood to be one-sided.) The points of discontinuity must satisfy: 

n + 1 + (7 — l)y 2 — i(2n + l)xy — 1(7 — l)xy 3 + ny 4 = 0. 

The last identity never holds for unitary x, y with n E N and < 7 < 2. For n = the right hand 
side of ([22]) is also continuous on Aj. 

2. For (x,y) E Aj x Ai we have (p(n,x,y) E Ai U A// U Aiv where An := {e i/3 |/5 E (7r/2,7r]} and 
A/y := {e^|/3 E [— 7r/2,0)} . To show this it is convenient to use the following substitution: 

u„ = tan — 



2 i{x n + 1) 
In the u-coordinates, (^) takes the form 

(u n + l)(n n _iPi(n,n n ) + P 2 («,it n )) 



«n+i = F(n,u n -i,u n ) := 



(u„ - l)(n„_iP 3 (n,n n ) + P±{n,u n )y 



where 



P 1 (n, v) = (2n + <y)v 3 - (2n + 4 + 7)1; 2 + (2n + 4 + 7)7; - (2n + 7), 
P 2 (n, u) = -(2n + 7)u 3 + (6n + 4 - j)v 2 + (6n + 4 - 7)1; - (2n + 7), 

p 3 ( n , w ) = (2n + 7)^ 3 + (6n + 4 - 7)w 2 - (6n + 4 - 7)?; - (2n + 7), 
P 4 (n, v) = -(2n + 7)u 3 - (2n + 4 + 7)v 2 - (2n + 4 + *y)v - (2n + 7). 
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Identity (|2^) reads as 

u = (n + l)(7ng-4-uo + 7 ) 

1 (n -l)(7«2 + 4no + 7)' 

We have to prove that for (u,v) £ [0,1] x [0,1], the values F(n,u,v) lie in the interval [— l,+oo]. 
The function F(n, u, v) is smooth on (0, 1) x (0, 1) and has no critical points in (0, 1) x (0, 1). Indeed, 
for critical points we have ^^(?v^, , '^ , ) = q which yields Pi(n,v)P^(n,v) — P2{n,v)P^{n,v) = and, 
after some calculations, v = 0,1,-1. On the other hand, one can easily check that the values of 
F(n, u, v) on the boundary of [0, 1] x [0, 1] lie in the interval [—1, +oo]. 

For n = 0, using ( |24| ) and exactly the same considerations as for F(n,0, v), one shows that — 1 < 
ui < +oo for uq £ [0, 1]. 

Now let us introduce 

Sn(k) := {x £ A^Xk £ A H ,X[ £ Ai V7 < I < k}, 

Siv(k) := {x £ Ai\x k £ A IV ,xi £ Ai V/ < I < k}, 

where x n is the solution of (|2C|). From the property 1 it follows that Sji(k) and Sjv(k) are open 

sets in the induced topology of Ai. 

Denote 

S n = uS n (k), Siv = uS IV (k), 

which are open too. These sets are nonempty since Su(l) and S7y(l) are nonempty. Finally 
introduce 

S/ := {x £ Aj : x n £ Ai \/n £ N}. 
It is obvious that Sj, Sji, and Siv are mutually disjoint. Property 2 implies 

5/U5/;U5/v = Ai. 

This is impossible for Sj = 0. Indeed, the connected set Aj cannot be covered by two open disjoint 
subsets Sn and Sjy. So there exists xo such that the solution x n £ Aj Vn. From 

(p(n, x, 1) = — i, (p(n, x, i) = — 1, (25) 

it follows that (for this solution) x n ^ 1, x n ^ i. This proves the theorem. 

To complete the proof of Theorem |l| it is necessary to show e* 77r//4 £ Si. This problem can be treated 
in terms of the method of isomonodromic deformations (see, for example, |7]] for a treatment of 
a similar problem). One could probably compute the asymptotics of solutions x n for n — > oo as 
functions of xq and show that the solution with x$ ^ e 77r / 4 cannot lie in Sj. The geometric origin 
of equation (]20| ) allows us to prove the result using just elementary geometric arguments. 

Proposition 4S / = {e^/ 4 }. 
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Proof: We have shown that Sj is not empty. Take a solution x n £ Si and consider the corresponding 
circle pattern (see Theorem |] and Theorem]^). Equations ( JlO| ) and (15) for N = M make it possible 
to find R(N + iN) and R(N + i(N + 1)) in a closed form. We now show that substituting the 
asymptotics of R(z) at these points into equation ( |TTj ) for M = N + 1, for immersed / ra ,m ; one 
necessarily gets i?(i) = tan 

Indeed, formula (|l9|) yields the following representation in terms of the T- function: 

R[N + iN) = c{ 1 Y r(iY + ^ /2) 



r(iV + i- 7 /2) : 

where 



7 r(i- 7 /2) r9R , 
c(7) = 2r(i + 7 /2) - (26) 



From the Stirling formula |T| 



I2tt ( s 

one obtains 



12s Vs 



rW = A/-(i) + — + ) (27) 



/?(A^ + /A-) = r(,)iVT- 1 (l + 0(I)). (28! 



Now let i?(i) = a tan ^ where a is a positive constant. Equation ( |l5|) for M = N, N > reads 

i?(iV - 1 + iN)R(N + i(iV + 1)) = i? 2 (iV + »JV). 

This is equivalent to the fact that the centers of all the circles C{N + iN) lie on a straight line. 
This equation yields 

R(N + i(N + 1)) = fatan ^V""" f - 1) + 7K2(^ - 3) + 7 )(2(iV - 5) + 7 ). 

// y 4/ I, (2JV-7)(2(JV-2)- 7 )(2(Ar-4)- 7 )... ) 

Using the product representation for tanx, 

t an x - sinx X ( 1 ~ ^) "' ( 1 ~ I^O 7 ) 

cosx ( 1 4x 2 \ / 1 4x 2 \ /-i 4z 2 \ 

one arrives at 

R(N + i(N + 1)) = a(- 1 ) iY c( 7 )^- 1 ) (l + o(^y (29) 
Solving equation (|Tl|) with respect to R 2 (z) we get 

R 2 {z) = G(N, M, R(z + i),R(z + 1), R(z - i)) := 



R{z + i)R{z + l)R(z -i) + R 2 {z + 1) (^fR(z -i) + ^fR(z + i] 
R(z + 1) + ^fR(z + l ) + ^R(z~^i) 



(30) 
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For z E V, R(z + i) > 0, + 1) > 0, R(z — i) > 0, the function G is monotonia 

<9G 

> °> ^7 TT > 0> ^777 T > 0- 



+ *) - ' dR{z + 1) - ' dR(z - i 
Thus, any positive solution R(z), z£Vof ([[]) must satisfy 

R 2 {z) > G(N, M, 0, R(z + 1), R(z - i)). 



Substituting the asymptotics of R ( |28|) and (|29| ) into this inequality and taking the limit K — > oo, 
for iV = 2K, we get a 2 > 1. Similarly, for N = 2K + 1 one obtains -K > 1, and finally o = 1. 
This completes the proof of the Proposition and the proof of Theorem p]. 



Remark. Taking further terms from the Stirling formula (j27p, one gets the asymptotics for Z 7 

^_m(2^y(i +0 (jfiy (31) 

having a proper smooth limit. Here the constant 0(7) is given by (p^). 



Due to representation (]?]) the discrete conformal map Z 7 can be studied by the isomonodromic 
deformation method. In particular applying a technique of M one can probably prove the following 



Conjecture The discrete conformal map Z" 1 has the following asymptotic behavior 



7 V 2 / V V Vn 2 + m 2 

For < 7 < 2 this would imply the asymptotic embeddedness of Z 7 at n, m — > 00 and, combined 
with Theorem [j], the embeddedness^ of Z 7 : — * C conjectured in 0,0]. 

5 The discrete maps Z 2 and Log. Duality 

Definition [| was given for < 7 < 2. For 7 < or 7 > 2, the radius R(l + i) = 7/(2 — 7) of 
the corresponding circle patterns becomes negative and some elementary quadrilaterals around /o,o 
intersect. But for 7 = 2, one can renormalize the initial values of / so that the corresponding map 
remains an immersion. Let us consider Z 7 , with < 7 < 2, and make the following renormalization 
for the corresponding radii: R — ► ^-^-R. Then as 7 — > 2 — from below we have 

R(0) = ?—t^+0, R(l+i) = l, R(i) = ?—ttan^ -> -. 

7 7 4 7T 



2 A discrete conformal map f n ,m is called an embedding if inner parts of different elementary quadrilaterals 

(fn,m, fn + l,m, fn + l,m + l, /n,m+l) do not intersect. 



16 



Definition 4 Z 2 : Z 2 — ► R 2 = C is the solution of (Ojj wrai/i 7 = 2 an<i i/ie initial conditions 

2 

Z 2 (0,0) = Z 2 (1,0) = Z 2 (0, 1) = 0, Z 2 (2,0) = l, Z 2 (0,2) = -l, Z 2 (l,l)=i-. 



Figure 6: Discrete Z A 



In this definition, equations (||) are @ are understood to be regularized through multiplication 
by their denominators. Note that for the radii on the border one has R(N + iN) = N. 
Equation (16) has the symmetry R — > 



Proposition 5 Let R(z) be a solution of the system (It, 11) for some 7. Then R(z) 
solution of (1C, \T1\) with 7 = 2 



B{z) 



is a 



7- 



This proposition reflects the fact that for any discrete conformal map / there is dual discrete 
conformal map f* defined by (see 



+l,m 



f* 



1 



n+l,m 



fn 



n. 



m+1 



n.r 



1 



fn,i 



m+1 



fn 



(32) 



Obviously this transformation preserves the kite form of elementary quadrilaterals and therefore is 



(f 



well-defined for Schramm's circle patterns. The smooth limit of the duality (32) is 

1 

~~~7 r 

The dual of f(z) = z 2 is, up to a constant, f*(z) = log z. Motivated by this observation, we define 
the discrete logarithm as the discrete map dual to Z 2 , i.e. the map corresponding to the circle 
pattern with radii 

-RLog(^) = l/Rz*{z), 
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where i?^ 2 are the radii of the circles for Z 2 . Here one has i?Log(0) = oo, i.e. the corresponding 
circle is a straight line. The corresponding constraint (0) can be also derived as a limit. Indeed, 
consider the map g = 2 ^-Z'~ 1 — ^p-. This map satisfies (Q) and the constraint 



/ , 2 7\ o (ffn+l.m 9n,m)i.9n.m 9n—l,m) , ~ (ffn.m+l 9n,m)\9n,ni 9n,m—X) 
7 SWm H = 2n t r h 2m- 



7 7 (ffn+l ,m 9n—l,m) \9n,m+l 9n,m—l) 

Keeping in mind the limit procedure used do determine Z 2 , it is natural to define the discrete 
analogue of log z as the limit of g as 7 — > +0. The corresponding constraint becomes 

-1 \9n+X,m 9n,m)\9n,m 9n—l,m) . \9n,m+l 9n,m){9n,rn 9n,m— l) /oo\ 
l = n t r \-m r . (33) 

\9n-\-l,m 9n—\,m) \9n,m+l 9n,m—\) 



Figure 7: Discrete Log. 



Definition 5 Log is the map Log : Z 2 — > R 2 = C satisfying ^) and fl#3[ ) mi/i i/ie initial condi- 
tions 

Log(0, 0) = 00, Log(l, 0) = 0, Log(0, 1) = m, 
Log(2, 0) = 1, Log(0, 2) = 1 + m, Log(l, 1) = i|. 



The circle patterns corresponding to the discrete conformal mappings Z 2 and Log were conjectured 



by O. Schramm and R. Kenyon (see [17]), but it was not proved that they are immersed. 
Proposition 6 Discrete conformal maps Z 2 and Log are immersions. 
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Proof: Consider the discrete conformal map — — -^-Z 7 with <^ 7 <C 2. The corresponding solution x n 
of (p0|) is a continuous function of 7. So there is a limit as 7 — > 2 — 0, of this solution with x n G Aj, 
xo = i, and x\ = e ^ ne solution x n of ( |20| ) with the property x n G Aj satisfies x n ^ 1, 

x n 7^ i for n > (see (p5|)). Now, reasoning as in the proof of Proposition |3], we get that Z 2 is 
an immersion. The only difference is that R(0) = 0. The circle C(0) lies on the border of V, so 
Schramm's result (see [}16[] ) claiming that corresponding circle pattern is immersed is true. Log 
corresponds to the dual circle pattern, with R^ og (z) = l/R Z 2(z), which implies that Log is also an 
immersion. 

6 Discrete maps Z 1 with 7 ^ [0, 2] 

Starting with Z 7 , 7 G [0, 2], defined in the previous sections, one can easily define Z 7 for arbitrary 7 
by applying some simple transformations of discrete conformal maps and Schramm's circle patterns. 
Denote by 5 7 the Schramm's circle pattern associated to Z 7 , 7 G (0,2]. Applying the inversion 
of the complex plane z \— > t(z) = 1/z to S" 7 one obtains a circle pattern tS^, which is also of 
Schramm's type. It is natural to define the discrete conformal map Z -7 , 7 G (0,2], through the 
centers and intersection points of circles of tSj. On the other hand, constructing the dual Schramm 
circle pattern (see Proposition ||) for Z -7 we arrive at a natural definition of Z 2+7 . Intertwining 
the inversion and the dualization described above, one constructs circle patterns corresponding to 
Z 7 for any 7. To define immersed Z 7 one should discard some points near (n, m) = (0, 0) from the 
definition domain. 

To give a precise description of the corresponding discrete conformal maps in terms of the 
constraint (^) and initial data for arbitrarily large 7 a more detailed consideration is required. 
To any Schramm circle pattern S there corresponds a one complex parameter family of discrete 
conformal maps described in Q. Take an arbitrary point Pqo G C U 00. Reflect it through all 
the circles of S. The resulting extended lattice is a discrete conformal map and is called a central 
extension of S. As a special case, choosing = 00, one obtains the centers of the circles, and 
thus, the discrete conformal map considered in Section ||. 

Composing the discrete map Z 7 : Z 2 + — > C with the inversion t(z) = 1/z of the complex plane 
one obtains the discrete conformal map G(n,m) = T(Z"'(n,m)) satisfying the constraint (§) with 
the parameter 7^ = —7. This map is the central extension of r5 7 corresponding to Poo = 0. Let us 
define Z -7 as the central extension of rS-y corresponding to Poo = 00, i.e. the extension described 
in Section H The map Z~ 7 defined in this way also satisfies the constraint (|2|) due to the following 

Lemma 3 Let S be a Schramm's circle pattern and f°° : — > C and f° : Z;j_ — > C be its two 
central extensions corresponding to Poo = 00 and Poo = 0, respectively. Then f°° satisfies (W if 
and only if f° satisfies ^). 

Proof: If /°° (or / ) satisfies @, then f£° (respectively f® ) lie on a straight line, and so do 
fo° m (respectively /o m )- A straightforward computation shows that /^° and /° satisfy @ simul- 
taneously, and the same statement holds for fQ° m and /o m - Since d) is compatible with @ f° 
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(respectively f°°) satisfy (|2j) for any n, m > 0. 

Let us now describe Z K for K G N as special solutions of (Jl], |2[) . 

Definition 6 Z K : Z 2 -> R 2 = C ; w/iere if G N, is i/ie solution of (@, ||j m'i/t 7 = if and f/ie 
initial conditions 

Z K (n,m) = for ra + m < if - 1, (n,m)GZ 2 _. (34) 
Z x (if,0) = l, (35) 



Figure 8: Discrete Z^. 
The initial condition ( |34|) corresponds to the identity 

^£ r ( 2; = o) = o, fc<if, 

in the smooth case. For odd K = 2N + 1, condition ( |36|) reads 



*W+l/o*nN_..(2tf-l)" 



Z /iV+i (2iV,l) 



(2iV)!! ' 



and follows from constraint (||). For even K = 2N, any value of Z K (K — 1,1) is compatible with 
(|2|). In this case formula ( |36|) can be derived from the asymptotics 

R(N + iN) 
hm — 



tv^oo R(N + i(N + 1)) 
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and reads 



r2N ln ,r , , ,2(2JV-2)!! 



Z 2iy (2N -1,1) =i— 



vr (2N - 1)!!' 

We conjecture that so denned Z K are immersed. 

Note that for odd integer K = 2iV+l, discrete Z 2N+1 in Definition || is slightly different from the 
one previously discussed in this section. Indeed, by intertwining the dualization and the inversion 
(as described above) one can define two different versions of Z 2N+1 . One is obtained from the circle 
pattern corresponding to discrete Z(n, m) = n + im with centers in n + im, n + m = (mod 2). The 
second one comes from Definition || and is obtained by the same procedure from Z{n, m) = n + im, 
but in this case the centers of the circles of the pattern are chosen in n + im, n + m = 1 (mod 2). 
These two versions of Z 3 are presented in Figure ||. The left figure shows Z 3 obtaied through 
Definition ||. Note that this map is immersed, in contrast to the right lattice of Figure ||] which has 
overlapping quadrilaterals at the origin (see Figure ||) . 




Figure 9: Detail view of two versions of discrete Z 3 . 
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Appendix A 

Proof of Theorem |^. 

Compatibility. Suppose that j n ,m is a solution of (|T|) which satisfies (H) for (n, 0) and (0, m) Vn, m. 
Direct, but very long, computation (authors used Mathematica computer algebra to perform it) 
shows that if the constraint (^) holds for 3 vertices of an elementary quadrilateral it holds for the 
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fourth vertex. Inductive reasoning yields that f n>m satisfies (|6|) for any n, m. So the constraint @ 
is compatible with (|l|). 

Necessity. Now let / n , m be a solution to the system (|l|),©. Define ^0,0 W as a nontrivial solution 
of linear equation (^) with A (A) given by Theorem ||. Equations (Q) determine ^ , n , m (A) for any 
re, rre. By direct computation, one can check that the compatibility conditions of (0) and @) 



d 



dX 
d 

dX 



Vr, 



(37) 



are equivalent to (jy, ^). 
Sufficiency. Conversely, let VPn.m 

(A) satisfy (0) and (|) with some A — ind.epGnd.ent matrices B nrn: 
C n ,m, -D n ,m- From (||) it follows that ti B 

n.m — i-i tr C n ,m — Equations ( |37| ) are equivalent 
to equations for their principal parts at A = 0, A = — 1, A = l, A = oo: 



D 



n+l,m 



D 



n,m+l 



1 



1 





B 



n+l.m I 
Bn,m+1 
Cn+l,m I 

c ( 



1 



1 



t'n 



u n,m 
1 



1 
1 



1 
1 



(-^n+l,m ^n+l,m Cn+l,m 
(D n ,rn+1 -^ri/m+l C nm _|_i 

From ([40], 0) and trS n . 77i — — ti, it follows that 








^ 





1 


:)- 


1 1 










^ 





1 


:)- 


I 1 






1 
1 



Dr. 



Dr. 



-M r , 



D r , 



1 



1 



1 



en 



1 



(D 
(D 



ft. TO B nm 



11. m B n m 



Cn,mt 

\c 

Cn,m) 
Cn,rn) 





1 



1 



(38) 
(39) 
(40) 
(41) 
(42) 
(43) 

(44) 
(45) 



n — <f> 



V"n,m U n m U n —\ m \ <f> 



U n ,m 4" Un- l.m \ 1 u n—l,m 



I. 
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Similarly, (fi^ , |4^) and tr C rhm = —m imply 

_ _ m-ip I v njm V„, m fn,-m-l \ _i^ T 
^n,m — I -i Jo 

v n,m H" v n,m—l V 1 "n,m- 1 / ^ 

Here, and ■0 are constants independent of n, m. The function a(A) in (|8|), independent of n and 
m, can be normalized to vanish identically, i.e. trD nm = 0. Substitution of 



D 

into equations (38, ^9|) yields 



a b 
c —a 



0"7i+l,m — Q"n,m Un,mCn,mi / n,m+l — Q"n,m ^n,m^n,m^ (^"0 
^n+l,m — ^n,m "I - 'U"n,m{.Q"n,m ~\~ Q"n+l.m)i &n,m+l — ^ra,m "I - ^ra,m(^ra,m "I - Q-n,m+l)- (^^) 

Thus c is a constant independent of re, m. Equations (fiT]) can be easily integrated 

where 6 is independent of n, m (recall that n n , m = f n +i,m-fn,m, «n,m = f n ,m+i-fn,m)- Substituting 
this expression into (E3) and integrating we get 



J n,m 



c fn,m 29fn,m + 



for some constant {l. Now (E3) and (45) imply 



n — (j) m — ip 

11"n,m ~r ^-n— l,m ^n,m \ ^n,m— 1 

which is equivalent to the constraint (||) after identifying c = § , = — J, fi = — § . 

Appendix B 

The proof of Lemma [2] uses the following technical 
Lemma 4 For positive R, the following hold: 



1) equations (|71y and ai z and equation (It) at z — i imply pqj at z + 1, 



2) equation ( \ldi ) at N + ziV and equations ( ]T(\ ) at N + iiV and ai iV — 1 + iiV imply 
(N + M){R(z) 2 - R(z + l)R(z - i)){R{z - i) + #0 - 1))+ 
(TV - M){R(zf - R{z - i)R(z - l))(R(z + 1) + i2(z - »)) = 0, (49) 
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at z = N + iM, for M = N + 1 
3) equations ( 



4) equations ( 



) and (4§) at z = N + iM, N / ±M, imply fcjjy at z, 
) and 0) at z = N + iM, N / ±M ; impZy 

(N + M)(R{zf - R(z + - 1))(.R(* + 1) + R(z + »))+ 

(iV - M)(R{zf - R{z + + i)){R{z + i) + - 1)) = 0, 



(50) 



and (11) at z, 



5) equations ( \5(\ ) and \l(\ ) at z and equation ( \l(\ ) at z — 1 imply ( mQ ) at z + i. 
Proof is a direct computation. Let us check, for example, 3). Equations (|l5|) and ( fi"9| ) read 
i{R{zf - R{z + i)R{z - 1)) + r](R(z + i) + R(z - 1)) = 0, 

i{R{zf - R{z + l)R(z - i)) - r](R(z + 1) + R(z - i)) = 0, 

where f = (iV + M)(22(z - 1) + - i)), rj = (M — N)(R{z) 2 - R{z - l)R(z - %)). Since f ^ we 
get 

(i?(z) 2 - + - l)){R{z + 1) + - i))+ 

(R(z) 2 - R{z + l)R(z - i))(R(z + i) + R{z - 1)) = 0, 
which is equivalent to fli~6|). 

Proof of Lemma [2[ 

By symmetry reasons it is enough to prove the Lemma for N > 0. Let us prove it by induction on 
N. 

For N = 0, identity @ yields 1 + iM) = + iM). Equation ^) at z = iM implies (||). 
Now equations (11) and (15) at z = iM imply (|l6|). 



Induction step N —> N + 1. The claim 1) of Lemma || implies ( p^5[ ) at z = N + 1 + iM. The claims 
2),3),4) yield equations (||) and (JTl|) at z = N + 1 + i(N + 2). Now using 5),3),4) of Lemma | one 
gets, by induction on L, equations (|T6|) and (O) at z = N + 1 + z(./V + L + 1) for any L £ N. 
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